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Abstract 


The strength and tensor structure of the Higgs boson’s interactions are investigated using an 
effective Lagrangian, which introduces additional CP-even and CP-odd interactions that lead to 
changes in the kinematic properties of the Higgs boson and associated jet spectra with respect 
to the Standard Model. The parameters of the effective Lagrangian are probed using a tit to five 
differential cross sections previously measured by the ATLAS experiment in the H —» yy decay 
channel with an integrated luminosity of 20.3 fb -1 at yfs = 8 TeV. In order to perform a simul¬ 
taneous fit to the five distributions, the statistical correlations between them are determined by 
re-analysing the H —> yy candidate events in the proton-proton collision data. No significant 
deviations from the Standard Model predictions arc observed and limits on the effective Lag¬ 
rangian parameters are derived. The statistical correlations are made publicly available to allow 
for future analysis of theories with non-Standard Model interactions. 
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1 Introduction 


The discovery of a Higgs boson at the ATLAS and CMS experiments [1,2] offers a new opportunity to search 
for physics beyond the Standard Model (SM) by examining the strength and structure of the Higgs boson’s 
interactions with other particles. Thus far, the interactions of the Higgs boson have been probed using the 
/r-framework [3], in which the strength of a given coupling is allowed to vary from the SM prediction by a 
constant value. In this approach, the total rate of a given production and decay channel can differ from the 
SM prediction, but the kinematic properties of the Higgs boson in each decay channel are unchanged. 

An alternative framework for probing physics beyond the SM is the effective field theory (EFT) approach [3- 
8], whereby the SM Lagrangian is augmented by additional operators of dimension six or higher. Some of 
these operators produce new tensor structures for the interactions between the Higgs boson and the SM 
particles, which can modify the shapes of the Higgs boson kinematic distributions as well as the associated 
jet spectra. The new interactions arise as the low-energy manifestation of new physics that exists at energy 
scales much larger than the partonic centre-of-mass energies being probed. 

In this Letter, the effects of operators that produce anomalous CP-even and CP-odd interactions between 
the Higgs boson and photons, gluons, W bosons and Z bosons are studied using an EFT-inspired effective 
Lagrangian. The analysis is performed using a simultaneous fit to five detector-corrected differential cross 
sections in the H —> yy decay channel, which were previously published by the ATLAS Collaboration [9]. 
These are the differential cross sections as functions of the diphoton transverse momentum (Pj 7 ), the number 
of jets produced in association with the diphoton system (Nj ets ), the leading-jet transverse momentum (pj ), 
and the invariant mass (mjj) and difference in azimuthal angle (A (pjj) of the leading and sub-leading jets 
in events containing two or more jets. The inclusion of differential information significantly improves the 
sensitivity to operators that modify the Higgs boson’s interactions with W and Z bosons. To perform a 
simultaneous analysis of these distributions, the statistical correlations between bins of different distributions 
need to be included in the fit procedure. These correlations are evaluated by analysing the H —> yy candidate 
events in the data, and are published as part of this Letter to allow future studies of new physics that produces 
non-SM kinematic distributions for H —> yy. 

2 Higgs effective Lagrangian 

The effective Lagrangian used in this analysis is presented in Ref. [8]. In this model, the SM Lagrangian is 
augmented with the dimension six CP-even operators of the Strongly Interacting Light Higgs formulation [6] 
and corresponding CP-odd operators. The H —> yy differential cross sections are mainly sensitive to the 
operators that affect the Higgs boson’s interactions with gauge bosons and the relevant terms in the effective 
Lagrangian can be specified by 


-Cell = CyOy + CgOg + ChwOhW + C hbO HB 

+ CyOy + Cgdg + CHWOHW + CHbOhB, 

where c; and c; are ‘Wilson coefficients’ specifying the strength of the new CP-even and CP-odd interactions, 
respectively, and the dimension-six operators (9, are those described in Refs. [8, 10]. In the SM, all of the 
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Wilson coefficients are equal to zero. The O y and O y operators introduce new interactions between the Higgs 
boson and two photons. The O CJ and O CJ operators introduce new interactions between the Higgs boson and 
two gluons and the analysis presented in this Letter is sensitive to these operators through the gluon fusion 
production mechanism. The Ohw and Ohw operators introduce new HWW , HZZ and HZy interactions. 
The HZZ and HZy interactions are also impacted by Ohb and Ohb and, to a lesser extent, O y and O y . The 
analysis presented in this Letter is sensitive to the Ohw, Ohw, Ohb and Ohb operators through vector-boson 
fusion and associated production. 

Other operators in the full effective Lagrangian of Ref. [8] can also modify Higgs boson interactions. Com¬ 
binations of some of the CP-even operators have been constrained using global fits to experimental data from 
LEPandthe LHC [8,11,12]. 

3 Statistical correlations between differential distributions 

ATLAS [13] is a multipurpose particle physics detector with cylindrical geometry and nearly 47r coverage 
in solid angle. 1 The analysis is performed using proton-proton collision data at a centre-of-mass energy 
yfs = 8 TeV and an integrated luminosity of 20.3 ftr 1 . 

The object and event selections used to define the differential distributions are described in detail in Ref. [9]. 
The statistical correlations between the measured cross sections as a function of different distributions are ob¬ 
tained using a random sampling with replacement method on the detector-level data. This procedure is often 
referred to as ‘bootstrapping’ [14]. Bootstrapped event samples are constructed from the data by assigning 
each event a weight pulled from a Poisson distribution with unit mean. The five differential distributions are 
then reconstructed using the weighted events, and the signal yields in each bin of a differential distribution 
are determined using an unbinned maximum-likelihood fit of the diphoton invariant mass spectrum (full de¬ 
tails of the fit can be found in Ref. [9]). The procedure is repeated 10000 times with statistically independent 
weights and the correlation between two bins of different distributions is determined from the scatter graph 
of the corresponding extracted cross sections. The observed correlations between bins of the measured p y 
and Aj e ts cross sections are shown in Fig. 1 . 

The statistical uncertainties on the correlation due to the finite number of bootstrap samples ranges from 
0.5% to 1%. The statistical uncertainty on the correlations due to the finite number of events in data is 
determined to be less than 2% using the statistical overlap and variance of signal and background events in a 
mass window around the Higgs boson mass. In order to validate this approach, a set of pseudo-experiments 
was created from input conditions (with known correlations) chosen to be similar to those in data in terms 
of purity, kinematics and sample size. For each pseudo-experiment, a value for the correlation is determined 
using 10000 bootstrapped samples and compared to the input correlations. No bias due to the bootstrapping 
is observed in the central value obtained from 500 pseudo-experiments. 


1 ATLAS uses a light-handed coordinate system with its origin at the nominal interaction point (IP) at the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis points upward. 
Cylindrical coordinates (r, (f>) are used in the transverse plane, (j> being the azimuthal angle around the beam pipe. 
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Figure 1: Statistical correlations between the measured cross sections in bins of the diphoton transverse momentum 
and jet multiplicity distributions. The quoted uncertainties refer to the total statistical uncertainty due to the finite 
number of bootstrapped samples and the finite number of data events. 


As part of this Letter, the correlations computed above are made publicly available in HEPDATA [15], 
allowing the analysis to be repeated using alternative effective Lagrangians, complete EFT frameworks, or 
other models with non-SM Higgs boson interactions. 


4 Theoretical predictions 

The effective Lagrangian has been implemented in FeynRules [10]. 2 Parton-level event samples are produced 
for specific values of Wilson coefficients by interfacing the universal file output from FeynRules to the 
Madgraph5 [17] event generator. Higgs boson production via gluon fusion is produced with up to two 
additional partons in the final state using leading-order matrix elements. The 0-, 1- and 2-parton events 
are merged using the MLM matching scheme [18] and passed through the Pythia6 generator [19] to create 
the fully hadronic final state. Event samples containing a Higgs boson produced either in association with 
a vector boson or via vector-boson fusion are produced using leading-order matrix elements and passed 
through the Pythia6 generator. For each production mode, the Higgs boson mass is set to 125 GeV [20] and 
events are generated using the CTEQ6L1 parton distribution function and the AUET2 parameter set [21]. 
All other Higgs boson production modes are assumed to occur as predicted by the SM. 

Event samples are produced for different values of a given Wilson coefficient. The particle-level differential 
cross sections are produced using Rivet [22], The Professor method [23] is used to interpolate between 


2 The implementation in Ref. [10] involves a redefinition of the gauge boson propagators that results in unphysical amplitudes 
unless certain physical constants are also redefined. The original implementation did not include the redefinition of these physical 
constants. However, the impact of redefining the physical constants is found to be less than 1% on the predicted cross sections 
across the range of Wilson coefficients studied. The reiative change in the predicted Higgs boson cross sections as functions of 
the different Wiison coefficients is aiso found to agree with that predicted by the Higgs characterisation framework [16], with less 
than 2% variation across the parameter ranges studied. 
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these samples, for each bin of each distribution, and provides a parameterisation of the effective Lagrangian 
prediction. The parameterisation function is determined using 11 samples when studying a single Wilson 
coefficient, whereas 25 samples are used when studying two Wilson coefficients simultaneously. As the 
Wilson coefficients enter the effective Lagrangian in a linear fashion, second-order polynomials are used 
to predict the cross sections in each bin. The method was validated by comparing the differential cross 
sections obtained with the parameterisation function to the predictions obtained with dedicated event samples 
generated at the specific point in parameter space. 

The model implemented in FeynRules fixes the Higgs boson width to be that of the SM, T# = 4.07 MeV 
[3]. The cross sections are scaled by T#/(T// + AT), where AT is the change in partial width due to a specific 
choice of Wilson coefficient. The change in partial width is determined for each Higgs coupling using the 
partial-width calculator in Madgraph5 and normalised to reproduce the SM prediction from Hdecay [24]. 

The leading-order predictions obtained from Madgraph5 are reweighted to account for higher-order QCD 
and electroweak corrections to the SM process, assuming that these corrections factorise from the new phys¬ 
ics effects. The differential cross section as a function of variable X for a specific choice of Wilson coefficient, 
Ci is given by 

^ = V fcf tof 35 / fef G5 

dx Zf\dx) I dx ) n [ dx ) n=0 ’ 

where the summation j is over the different Higgs boson production mechanisms, ‘MG5’ labels the Mad- 
graph5 prediction and ‘ref’ labels a reference sample for SM Higgs boson production. 

The reference sample for Higgs boson production via gluon fusion is simulated using MG5_aMC@NLO [25] 
with the CT10 parton distribution function [26]. The H + //-jets topologies are generated using next-to- 
leading-order (NLO) matrix elements for each parton multiplicity (n = 0. I or 2) and combined using the 
FxFx merging scheme [27], The parton-level events are passed through Pythia8 [28] to produce the hadronic 
final state using the AU2 parameter set [29], The sample is normalised to the total cross section predicted 
by a next-to-next-to-leading-order plus next-to-next-to-leading-logarithm (NNLO+NNLL) QCD calculation 
with NLO electroweak corrections applied [3]. The reference sample for Higgs boson production via vector- 
boson fusion (VBF) is generated at NLO accuracy in QCD using the Powheg Box [30]. The events are 
generated using the CT10 parton distribution function (PDF) and Pythia8 with the AU2 parameter set. The 
VBF sample is normalised to an approximate-NNLO QCD cross section with NLO electroweak corrections 
applied [3]. The reference samples for Higgs boson production in association with a vector boson ( VH , 
V=W,Z) or a top-anti top pair (tJH) are produced at leading-order accuracy using Pythia8 with the CTEQ6L1 
PDF and the 4C parameter set [21]. The ZH and WH samples are normalised to cross sections calculated at 
NNLO in QCD with NLO electroweak corrections, whereas the ttH sample is normalised to a cross section 
calculated to NLO in QCD [3]. 

The ratio of the differential cross sections to the SM predictions for some representative values of the Wilson 
coefficients are shown in Fig. 2. The impact of the c g and c g coefficients are presented for the gluon fusion 
production channel and show a large change in the overall cross section normalisation. The c g coefficient 
also changes the shape of the A(pjj distribution, which is expected from consideration of the tensor structure 
of CP-even and CP-odd interactions [31,32]. The impact of the chw and chw coefficients are presented for 
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Figure 2: Ratio of differential cross sections predicted by specific choices of Wilson coefficient to the differential cross 
sections predicted by the SM. 


the VBF+VH production channel and show large shape changes in all of the studied distributions. 3 The 
A (pjj distribution is known to discriminate between CP-odd and CP-even interactions in the VBF production 
channel [34]. 


5 Limit-setting procedure 


Limits on the Wilson coefficients are set by constructing ax 2 function 


X ~ (cTdata f J ~p r ed j C ^ ( r l j a p l 0" pred ^ 


where (r daLa and <r pred are vectors from the measured and predicted cross sections of the five analysed ob¬ 
servables, and C = C st at + C exp + C pre d is the total covariance matrix defined by the sum of the statistical, 
experimental and theoretical covariances. The predicted cross section <T pred and its associated covariance 
C pre d are continuous functions of Wilson coefficients. Scans of one or two Wilson coefficients are earned 
out and the minimum x 2 value, x 2 ■ , is determined. The confidence level (CL) of each scan point can be 
calculated as 


1 - CL = 


X OO 

2 (ci)~xt 


d x f{x\ m ), 


3 Form factors are sometimes used to regularise the change of the cross section above a momentum scale A FF . This was investigated 
by reweighting the VBF+VH samples using form-factor predictions from VBFNLO [33]. The impact on the Chw and chw limits 
are negligible for A FF > 1 TeV. 
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Figure 3: The input data from Ref. [9] is compared to the SM hypothesis and two non-SM hypotheses with c g = 1x10 4 
and chw = 0.05, respectively. 


with^ 2 (c/j being the;y 2 value evaluated for a given Wilson coefficient c,, and /(x; m) being the;y 2 distribution 
for m degrees of freedoms and n = 1 or j for two-sided or one-sided limits. The coverage of CL and the 
effective number of degrees of freedom are determined using ensembles of pseudo-experiments. 4 

The input data vector is compared in Fig. 3 to the SM hypothesis as well as two non-SM hypotheses specified 
by c g ~ 1 x 1 0 -4 and chw = 0-05, respectively. 

The covariance matrix for experimental systematic uncertainties is constructed from all uncertainty sources 
provided by Ref. [9], which include the jet energy scale and resolution uncertainties, photon energy and 
resolution uncertainties, and model uncertainties. Identical sources are assumed to be fully correlated across 
bins and variables and the sign of an error amplitude is taken into account when computing the covariance 
matrix. The statistical uncertainties on the cross correlation have a negligible impact on the results reported 
here. 

The covariance matrix for the theoretical uncertainties is constructed to account for missing higher-order 
corrections and PDF uncertainties in the SM reference predictions. The uncertainties in the gluon fusion ref¬ 
erence samples are: (i) a shape uncertainty, estimated by simultaneously varying the factorisation and renor¬ 
malisation scales in MG5_aMC@NLO by a factor of 0.5 or 2.0, and (ii) uncertainties on the NNLO+NNLL 
QCD plus NLO electroweak total cross-section prediction [3], arising from missing higher-order corrections 
and PDF uncertainties; these uncertainties are assumed to be fully correlated among bins and observables. 

4 For one-dimensional limits on the CP-even (odd) Wilson coefficients, good agreement is found between the asymptotic formula 
and the pseudo-experiment test statistic with m = 1 and n = 1 (|). For the two-dimensional limits on c g versus c g , and Chw 
versus c alv , good agreement between pseudo-experiments and asymptotic formula is found for m = 1 and n = 1. For the two 
dimensional limit on c g versus c 7 , good agreement between pseudo-experiments and asymptotic formula is found for m = 2 and 

n = 1. 
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Figure 4: The 68% (dark) and 95% (light) confidence regions for the fit to the c r and c g Wilson coefficients. All other 
coefficients are set to zero. The shaded area represents the allowed region of parameter space and the marker indicates 
the SM value. 


For VBF, ZH and WH, shape uncertainties are neglected because their impact is expected to be negligible 
with respect to all other theory uncertainties. Normalisation uncertainties for these processes are taken from 
Ref. [3], 

The benefit of using more than one differential distribution in the analysis is quantified using an ‘Asimov 
dataset’, which is a representative dataset of the median expected cross-section measurement assuming the 
SM. For c g and c g , the use of a single inclusive distribution (/;: | ’ > or Ajets) results in the same expected limits 
as the full five-dimensional fit. For c y and c y , the most sensitive variable is found to be p y y , with a 5% 
improvement in the expected limits obtained from using the five-dimensional information. For chw an d 
chw , the most sensitive variable is A <pjj and an 18% improvement in the expected limits is obtained from 
using the five-dimensional fit. In summary, the expected sensitivity for c g , c g , c y and c y arises mainly from 
the normalisation of the different production mechanisms, and can be probed using the inclusive distributions 
that distinguish between the different processes, whereas the chw and chw coefficients benefit more from 
the full five-dimensional information due to the induced shape changes in the kinematics of the VBF+VF1 
process. 


6 Results 

The 68% and 95% confidence regions for a two-dimensional scan of c y and c g are shown in Fig. 4, after 
setting all other Wilson coefficients to zero. These additional interactions can interfere with the correspond- 
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Figure 5: The 68% (dark) and 95% (light) confidence regions for the fit to the c g and c CJ Wilson coefficients. All other 
coefficients are set to zero. The shaded area represents the allowed region of parameter space and the marker indicates 
the SM value. 


ing SM interactions. Destructive interference, for example, causes the H —> yy branching ratio to be zero at 

c y ~ 2 x 10 3 and the gluon fusion production cross section to be zero at c g -2.2 x 10 4 . The impact of 

these effects is evident in the structure of the obtained limits in the two-dimensional parameter plane. 

The 68% and 95% confidence regions for a two-dimensional scan of c g and c g are shown in Fig. 5, after 
setting all other Wilson coefficients to zero. The Ac fijj distribution is sensitive to the c g parameter through 
the gluon fusion production mechanism (Figs. 2 and 3) and the limit on c g is improved with the inclusion of 
this data in the fit. This is evident in Fig. 5 where the limit band is constricted at the largest values of c g . 

The 68% and 95% confidence regions obtained from scanning cuw and cuw are shown in Fig. 6, after 
setting cub — chw and c HB = cuw to ensure that the partial width for H Zy is unchanged from the 
SM prediction. 5 As discussed in Section 5, these Wilson coefficients produce large shape changes in all 
distributions and the obtained limits are strongest when fitting all five distributions simultaneously. 

The 95% confidence regions for cuw and cuw can be translated into the Higgs Characterisation frame¬ 
work [16] and compared to the ATLAS results for non-SM CP-even and CP-odd HVV interactions, which 
were obtained using an angular analysis of the decay products in the WW* and ZZ* decay channels [36]. 
The translated limits are -0.08 < £hvv/%sm < 0 09 and -0.22 < tanfcr) • Lvv/^sm < 0.22, where the 
variables /?hvv- £avv< Tsm and a are defined in Refs. [16,36]. The limits obtained in this analysis are a factor 
of approximately seven stronger than those in Ref. [36], due to increased sensitivity to the different Higgs 

5 Values of | chw ~ £hb\ > 0.033 lead to a very large decay rate for the H —» Zy process that is contradicted by the experimental 
constraints reported by ATLAS [35]. 
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C HW 


Figure 6: The 68% (dark) and 95% (light) confidence regions for the fit to the chw and chw Wilson coefficients. All 
other Wilson coefficients are set to zero, except for chb and chb which are set to be equal to chw and chw, respectively. 
The shaded area represents the allowed region of parameter space and the marker indicates the SM value. 


boson production channels arising from the inclusion of rate and jet kinematic information in the signal 
hypothesis. 

The observed limits on chw and chw arc also not excluded by current signal strength measurements. For 
example, the signal strength in the H —> ZZ* and H —> WW* channels is predicted to be approximately 1.3 
for chw - 0.1, which is consistent with the dedicated measurements [37,38]. 

The 95% confidence regions for a one-dimensional scan of the Wilson coefficients are given in Table 1. 


7 Summary 


The strength and structure of the Higgs boson’s interactions with other particles have been investigated us¬ 
ing an effective Lagrangian. Limits arc placed on anomalous CP-even and CP-odd interactions between 
the Higgs boson and photons, gluons, W -bosons and Z-bosons, using a fit to five differential cross sections 
previously measured by ATLAS in the H —» yy decay channel at -\[s = 8 TeV [9], No significant deviations 
from the SM predictions are observed. To allow a simultaneous fit to all distributions, the statistical correla¬ 
tions between these distributions have been determined by re-analysing the candidate H —» yy events in the 
proton-proton collision data. These correlations are made publicly [15] available to allow for future analysis 
of theories with non-SM Higgs boson interactions. 
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Coefficient 

95% 1 - CL limit 

Cy 

[-7.4, 5.7]xl0~ 4 U [3.8, 5.1]xl0” J 

Cy 

[-1.8, 1.8]xlO~ 3 

Cg 

[-0.7, 1.3]xl0' 4 U [-5.8, -3.8]xl0~ 4 

Cg 

[-2.4, 2.4]xl0~ 4 

chw 

[-8.6, 9.2]xl0~ 2 

chw 

[-0.23, 0.23] 


Table 1: Observed allowed ranges at 95% CL for the c r , c g and chw Wilson coefficients and their CP-conjugate partners. 
Limits on c y , c g , c y and c g are each derived with all other Wilson coefficients set to zero. Limits on chw and chw are 
derived with chb = chw and chb — chw , respectively. Two allowed regions are observed for c y and c g , with the region 
between the solutions producing too small pp —> H —> yy cross section due to destructive interference between new 
interactions and the SM. 
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